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ABSTRACT 

We study the BL Lac objects detected in the one year all sky survey of the Fermi satellite, with 
a energy spectral slope in the [0.1-100 GeV] band greater than 1.2. In the vs 7-ray 
luminosity plane, these BL Lacs occupy the region populated by Flat Spectrum Radio Quasars 
(FSRQs). Studying the properties of their spectral energy distributions (SED) and of their 
emitting lines, we find that several of these BL Lacs have a SED similar to FSRQs and that 
they do have broad lines of large equivalent width, and should be reclassified as FSRQs even 
adopting the current phenomenological definition (i.e. equivalent width EW of the emitting 
line greater than 5 A). In other cases, even if the EW width is small, the emitting lines can 
be as luminous as in quasars, and again their SED is similar to the SED of FSRQs. Sources 
classified as BL Lacs with a SED appearing as intermediate between BL Lacs and FSRQs also 
have relatively weak broad emission lines and small EW, and can be considered as transition 
sources. These properties are confirmed also by model fitting, that allows to derive the relevant 
intrinsic jet parameters and the jet power This study leads us to propose a physical distinction 
between the two classes of blazars, based on the luminosity of the broad line region measured 
in Eddington units. The dividing line is of the order of iBLR/^Edd ^ 5 x 10^'', in good 
agreement with the idea that the presence of strong emitting lines is related to a transition in 
the accretion regime, becoming radiatively inefficient below a disk luminosity of the order of 
one per cent of the Eddington one. 

Key words: BL Lacertae objects: general — quasars: general — radiation mechanisms: non- 
thermal — gamma-rays: theory — X-rays: general 



1 INTRODUCTION 

Among the blazars detected by the Large Area Telescope (LAT) 
onboard the Fermi satellite after 1 1 months of all sky survey (Abdo 
et al. 2010a, hereafter AlO) there are roughly an equal number 
of sources identified as BL Lac objects and Flat Spectrum Radio 
Quasars (FSRQs). The corresponding catalog of AGN detected at 
high Galactic latitude (|&| > 10°) is called First LAT AGN Cata- 
log (ILAC). In general, the LAT-detected BL Lac objects and Flat 
Spectrum Radio Quasars (FSRQs) separate quite well in the 7-ray 
spectral index - 7-ray luminosity plane (a-y — L^, where is the 
energy spectral index), in agreement with the early results borne out 
with the 3-months all sky survey of Ferm/ZLAT, which contained 
only 1/7 of the blazars in the ILAC catalog (Ghisellini, Maraschi 
& Tavecchio 2009, hereafter GMT09). On the other hand, there 
are a number of sources, classified as BL Lac objects, located in 
the region of the plane preferentially "inhabited" by FSRQs: these 
are BL Lacs with a relatively steep spectrum (i.e. a-y > 1.2). These 
"intruders" has been classified as BL Lac objects on the basis of the 
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"historical" distinction among BL Lacs and FSRQs, i.e. by means 
of the equivalent width (EW) of their emission lines (see e.g. Urry 
& Padovani 1995). Objects with a rest frame EW<5 A are called 
BL Lacs. This definition has the obvious advantage of being simple 
and of immediate use for an observational characterization of the 
object. On the other hand, the optical continuum of most blazars is 
relativistically enhanced by beaming, and very variable. In several 
cases a small EW does not imply emission lines of low luminos- 
ity, being simply the result of a particularly beamed non-thermal 
continuum. On the opposite side, EW greater than 5 A may be the 
results of a particularly low state of the beamed continuum in a 
source of intrinsically weak lines. A division based on the EW of 
emission lines does measure the relative importance of the beamed 
non-thermal continuum and the underlying thermal emission, but 
after the discovery that most of the non-thermal emission is at 7- 
ray energies, we know that the optical non-thermal flux very often 
is a minor contribution to the total, bolometric, non-thermal out- 
put. Therefore the EW alone is not a good indicator of the relative 
importance of the two contributions. 

Up to now, we construct samples of BL Lac as well as of 
FSRQs in order to study their properties and their possible differ- 
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ences, and adopt the classical, EW-based, sub-division. If the aim 
is to study intrinsically physical properties, this may be dangerous, 
since with the EW classification we may - for instance - classify 
as a BL Lac object a source with very luminous lines, typical of a 
FSRQ, only because at the time of the spectroscopic observations 
leading to the measurement of the EW the optical non-thermal flux 
was particularly intense. For illustration, let us take the case of 
PKS 0208-512. It has an observed Mgll emission line of EW~5 
A (2.5 in the rest frame), whose luminosity is close to 10** erg 
s~^, stronger than in some FSRQs. This object is classified as a BL 
Lac, but all its physical properties are resembling FSRQs. 

We therefore believe that a new classification scheme is 
needed, based on a physical property of the source. We suggest 
a division based on the luminosity of the broad emission lines, 
normalized to the corresponding Eddington luminosity, the natu- 
ral luminosity-scale. Normalizing in this way allows to compare 
objects of different black hole masses. This division implies to es- 
timate the black hole mass, that it is not a direct observable quantity. 
On the other hand, in recent years, the establishing of correlations 
between i) the luminosity of the bulge of the host galaxy and the 
black hole mass (Magorrian et al. 1998; Bentz et al. 2009); ii) the 
correlation between the dispersion velocity and the black hole mass 
(Ferrarese & Merritt 2000; Gultekin et al. 2009), and iii) the cor- 
relation between the luminosity of the continuum at selected fre- 
quencies and the size of the Broad Line Region (BLR, Vestergaard 
2002; Decarli et al. 2010 and references therein), made the esti- 
mate of the black hole mass much more affordable. Furthermore, 
in specific cases, very powerful blazars do have their IR-optical- 
UV continuum dominated by a thermal component produced by 
their accretion disk: modelling it with a standard Shakura-Sunyaev 
(1973) disk allows to find both the black hole mass and the accre- 
tion rate. 

We then investigate if the "intruder" BL Lacs in the a-y-L^ 
plane have intrinsically weak emission lines (in Eddington units) or 
if instead their EW is only a consequence of a particularly enhanced 
non-thermal continuum, or else if they are transition objects, with 
intermediate values of the broad line luminosity. 

We use a cosmology with h = SIa = 0.7 and SIm = 0.3, 
and use the notation Q — 10^ in cgs units (except for the black 
hole masses, measured in solar mass units). 



2 THE FERMI BLAZARS' DIVIDE ONE YEAR AFTER 

Fig. [U shows the energy spectral index a-y as a function of the (K- 
corrected, see GMT09) 7-ray luminosity for the FSRQs and 
BL Lacs on the I LAC sample of known redshift. This figure can 
be compared with the same one in GMT09 reporting the bright 
blazars of the LAT Bright AON Sample (LBAS sample (Abdo et al. 
2009, hereafter A09) for the 3-months all sky survey. In that figure 
there was a specific 7-ray luminosity dividing FSRQs and BL Lacs, 
around a few times 10*" erg s^^, interpreted as a consequence 
of the changing accretion regime of the underlying accretion disk 
from radiatively efficient to inefficient, or, in other words, from 
a standard Shakura Sl Sunyaev (1973) disk to an ADAF (Advec- 
tion Dominated Accretion Flow) one. As remarked in GMT09, the 
appearance of the dividing luminosity required that all the bright 
blazars in the LBAS sample have black holes of the same mass and 
similar beaming factors. This is approximately appropriate when 
considering the brightest sources, but when the decreased limiting 
flux allows to explore smaller luminosities both for FSRQs and 
BL Lacs, then it is likely that the corresponding black hole mass 



(and/or the beaming factor) is smaller, and the dividing luminos- 
ity gets spread into a larger range of values (as large as the spread 
in black hole masses and beaming factors). This explains why, in 
Fig.[T] no dividing luminosity is present. A comparison with Fig. 1 
of GMT09 shows that all the most luminous blazars were present 
also in LBAS: the decreased flux limit did not led to discover any 
new more powerful object. Furthermore, differently from Fig. 1 of 
GMT09, there is now no trend between the minimum luminosities 
and spectral index, although, at low luminosities, flat 7-ray spectra 
BL Lacs are more numerous than sources with a relatively steep a-y . 
This is likely due to LAT being more sensitive to flat spectra than 
to steep ones (see Fig. 9 of AIO). As already noted in GMT09, a 
correlation between and Ly is not expected, since at low 7-ray 
luminosities we expected the detection of FSRQs of lower black 
hole masses and smaller beaming. This fills the top left part of the 
a-y-L-y plane. 

The majority of BL Lac sources are characterized by a rela- 
tively flat (ay<1.2), but there are several exceptions. Some of 
these BL Lacs, however, have been classified as such because the 
equivalent width (EW) of their broad lines (that are indeed present) 
is less than 5 A. We discussed in GMT09 the cases of PKS 0537- 
441, AO 0235-1-164 and PKS 0426-380, that do have broad lines 
visible in their low emission states (see Sbarufatti et al. 2005 for 
PKS 0426-380; Plan et al. 2002 for PKS 0537^41; Raiteri et al. 
2007a for AO 0235-I-I64). Another example is 0208-512: it was 
observed to have a Mgll broad line with EW~5A (2.5A in the rest 
frame), but with a very large luminosity (~ 10** erg s^^; Scarpa 
& Falomo 1997). Therefore 0208-512 (and the other mentioned 
sources) are FSRQs whose non-thermal continuum is enhanced so 
much to make the very luminous broad lines to almost disappear, 
and not BL Lac objects with genuinely weak lines. Not appreciat- 
ing this point may cause some confusion when comparing FSRQs 
and BL Lacs. 

With respect to the LBAS sample, the number of sources clas- 
sified as BL Lacs but of unknown redshift increased: from the 
source list in AIO there are 159 sources classified as "BLL" with 
no redshift (in the "clean" sample), excluding the sources classified 
as of "unknown" type. 

In Fig. [2] we show the energy spectral index a-y vs the 0. 1- 
100 GeV photon flux of these 159 "BLL" sources. Of these 159 
sources, 54 have a-y ^ 1.2 (80 if we include the sources with "un- 
certain" classification). Fig.|2]shows in grey the area where most of 
the sources are. Exceptions at large photon fluxes are labelled, and 
we here briefly comment about these sources. 

PKS 1424-1-240 and PG I553-HI13 have been detected at TeV 
energies (see Ong et al. 2009; Teshima et al. 2009 for I424-I-240, 
and Aharonian et al. 2006; Albert et al. 2007 for PG 1553-1-1 13; see 
also Prandini et al. 2010), and very likely also 3C 66A (Acciari et 
al. 2009), although there can be a contamination from the nearby 
3C 66B radiogalaxy (see the discussion in Tavecchio & Ghisellini 
2009). The redshift of 3C 66 A is uncertain, even if a value of z = 
0.444 is commonly used. Due to the TeV detection, these 3 sources 
cannot lie at very large redshift although z up to ^--^0.7 would be 
possible (see Tavecchio et al. 2010 for PKS 1424-1-240 modelled 
assuming z = 0.67). 

For B3 0814-1-425 the NED database gives z = 0.53, quoting 
Sowards-Emmerd et al. (2005) from SDSS data. However, the in- 
spection of the SDSS spectrum does not confirm this redshift (nor 
the other quoted value, 2 = 1.07). 

Finally, there is no information for the redshift of CRATES 
J1542-I-6I29 = GB6 11542-^6129. It has been imaged by the SDSS, 
but no spectrum is available. 
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Table 1. in the LAT band (0.1-100 GeV) in units of 10 * ph cm s ^. L^, in the same band, is k-corrected and in units of erg s ^. *: no Swift 
observations. Sources whose name is in italics are present in Ghisellini et al. (2010a, hereafter GIO), and some of them are present in Tavecchio et al. (2010). 



If we assign to all sources in the grey area of Fig. [2] a given 
redshift, we can see the corresponding region in Fig. [T] We show 
this for two redshifts: z = 0.5 and z = 2, as labelled. It can be seen 
that in the case of z ~ 0.5 the BLL sources of unknown redshift 
would lie in the region already occupied by the other BL Lacs, 
while they would be "outliers" if the redshift is as large as 2 (see 
also the discussion in AlO). 

In other words: if the BLL sources in the ILAC catalog with 
unknown redshift will turn out to be at z<0.5-l, then they will fit 
in the phenomenological blazar sequence [i.e. they would be BL 
Lacs of low and moderate luminosity, with the majority having a 
flat 7~ray slope (i.e. < 1-1.2); Fossati et al. 1998; Donato et al. 
2001], while they would pose a problem if their redshift is larger. 
We will further discuss this point in §5). 



2.1 The "intruder" BL Lacs sample 

We consider all the sources classified as BL Lacs in the "clean" 
ILAC sample (defined as BL Lacs with \h\ > 10°, detected with 
a TS significance larger than 25 (TS stands for Test Statistics, see 
Mattox et al. 1996 for the definition. TS=25 approximately corre- 
sponds to 5(j), whose identification is secure and unique. We se- 
lected the sources of known redshift with a 7-ray energy spectral 
index a-y larger than 1.2, corresponding to a photon spectral index 
T-y > 2.2. The resulting 28 BL Lac objects are listed in Tab.[T] At 



the end of the same table we add the two BL Lacs (0235+164 and 
0426-380) that were "intruders" because of their extremely large 
7-ray luminosity (i.e. Lj > 10^** erg s^^), even if their spectral 
index was somewhat flatter than a-, — 1.2. All these objects are 
shown and labelled in Fig.[T] and correspond to the filled squares. 

We will characterize the SED of all the 30 BL Lacs of Tab.^ 
and search for existing data of the luminosity of their broad lines, 
if present. 



3 SWIFT OBSERVATIONS AND ANALYSIS 

Several blazars studied in this paper were observed by the Swift 
satellite. These objects are listed in Tab. 2 (X-ray data) and Tab. 
3 (optical/UV data). Even when they were performed during the 
1 1 months of the ILAC survey, they correspond to a "snapshot" of 
the optical-X-ray state of the source, while the 7-ray data are an 
average over the 1 1 months. Given the very rapid blazar variability, 
the SEDs constructed in this way should be considered, in all cases, 
not simultaneous (even when the Swift UVOT and XRT data are 
indeed simultaneous). 

The data were screened, cleaned and analysed with the soft- 
ware package HEASOFT v. 6.8, with the calibration database up- 
dated to 30 December 2009. The XRT data were processed with the 
standard procedures (XRTPIPELINE v. . 12 . 4). All sources 
were observed in photon counting (PC) mode and grade 0-12 (sin- 
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[erg s-i] 

Figure 1. The energy spectral index as a function of tiie 7-ray luminosity L-^ in the band [0. 1-10 GeV] for all blazars with known redshift present in the 
ILAC sample. The filled circles (red in the electronic version) are FSRQs; empty (blue) squares are BL Lacs with o-, < 1.2, and fileld (green) squares are 
sources classified as BL Lacs in the ILAC sample with a-, > 1.2. In addition, the two larger (cyan) circles are 0235+164 and 0426-380, classified as BL 
Lacs, that have > lO'** erg The horizontal grey line marks = 1.2. The two grey regions illustrate how the corresponding grey area shown in Fig. 
|2]would lie assuming a redshift of 0.5 or 2, as indicated. 



gle to quadruple pixel) were selected. The channels with energies 
below 0.2 keV and above 10 keV were excluded from the fit and 
the spectra were rebinned in energy so to have at least 20-30 counts 
per bin in order to apply the test. When there are no sufficient 
counts, then we applied the likelihood statistic in the form reported 
by Cash (1979). Each spectrum was analysed through XSPEC v. 
12.5. In with an absorbed power law model with a fixed Galactic 
column density as measured by Kalberla et al. (2005). The com- 
puted errors represent the 90% confidence interval on the spectral 
parameters. Tab. 2 reports the log of the observations and the best 



fit results of the X-ray data with a simple power law model. The 
X-ray spectra displayed in the SED have been properly rebinned to 
ensure the best visualization. 

UVOT (Roming et al. 2005) source counts were extracted 
from a circular region 5"-sized centred on the source position, 
while the background was extracted from a larger circular nearby 
source-free region. Data were integrated with the uvotimsum 
task and then analysed by using the uvot source task. The ob- 
served magnitudes have been dereddened according to the formulae 
by Cardelli et al. (1989) and converted into fluxes by using standard 
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Figure 2. Energy index a-y vs 7-ray flux in the [0.1-100 GeV] energy 
band for the sources classified in the clean ILAC catalog as "BLL" (and 
excluding the ones classified as "unknown"). The plotted [0.1-100 GeV] 
photon flux has been calculated from the [1-100 GeV] flux listed in AlO, 
using the corresponding spectral index. The grey region corresponds to the 2 
grey areas in Fig.[T] The (green) squares are the 56 sources with a-y ^ 1.2 
(this limit is shown by the horizontal line). See text for a brief comment 
about the labelled sources that are located outside the grey region. 



foimulae and zero points from Poole et al. (2008). Tab. |3] lists the 
observed magnitudes. 



4 MODELLING THE SED 

To model the SEDs of the blazars in this sample we used the same 
model used in Ghisellini et al. (2010a, hereafter GIO). It is a one- 
zone, leptonic model, fully discussed in Ghisellini & Tavecchio 
(2009). In that paper we emphasize the relative importance of the 
different sources of the seed photons for the inverse Compton scat- 
tering process, and how they change as a function of the distance of 
the emitting region from the blacli hole. Here we briefly summarize 
the main characteristics of the model. 

The source is assumed spherical (radius R) and located at a 
distance i?diss from the central black hole. The emitting electrons 
are injected at a rate (3(7) [cm~^ s"'^] for a finite time equal to the 
light crossing time R/c. The shape of (5(7) we adopt is assumed 
to be a smoothly broken power law with a break at 7b: 



Q(7) = Qo ■ 



{ihh 



(1) 



1 + (7/7b)-"i+=2 
The emitting region is moving with a velocity j3c corresponding 
to a bulk Lorentz factor T. We observe the source at the viewing 
angle 6^ and the Doppler factor is(5=l/[r(l— (5 cos 6^)]. The 
magnetic field B is tangled and uniform throughout the emitting 
region. For the inverse Compton process, besides the synchrotron 
seed photons (produced internally to the jet), we take into account 
several sources of radiation produced externally to the jet: i) the 
broad line photons, assumed to re-emit 10% of the accretion lumi- 
nosity from a shell-like distribution of clouds located at a distance 
i?BLR = 10^^1/^45 cm; ii) the IR emission from a dusty torus. 



located at a distance Rm. = 2.5 x IQ^^ L]/^^ cm and reprocessing 
10%-30% of the accretion luminosity; iii) the direct emission from 
the accretion disk, including its X-ray corona. We also consider 
the starlight contribution from the inner region of the host galaxy 
and the cosmic background radiation, but these photon sources are 
unimportant in our case. All these contributions are evaluated in 
the blob comoving frame, where we calculate the corresponding 
inverse Compton radiation from all these components, and then 
transform into the observer frame. 

We calculate the energy distribution (7) [cm^^] of the emit- 
ting particles at the particular time R/c, when the injection process 
ends. Our numerical code solves the continuity equation which in- 
cludes injection, radiative cooling and pair production and re- 
processing. Ours is not a time dependent code: we give a "snap- 
shot" of the predicted SED at the time R/c, when the particle distri- 
bution A'^(7) and consequently the produced flux are at their maxi- 
mum. 

The accretion disk component is calculated assuming a stan- 
dard optically thick geometrically thin Shakura & Sunjaev (1973) 
disk. The emission is locally a black body. The temperature pro- 
file of the disk is given e.g. in Frank, King & Raine (2002). In our 
sources, the optical-UV continuum is almost always dominated by 
the beamed non-thermal emission. On the other hand, when de- 
tected, the broad emission lines allow to estimate the luminosity of 
the accretion disk Ld- In these cases we have assumed the value of 
Ld derived from the emission lines. 

By estimating the physical parameters of the source we can 
calculate the power that the jet carries in the form of radiation (Pr), 
magnetic field (Pb), relativistic electrons (Pe) and cold protons 
(Pp) assuming one proton per electron. These powers are calcu- 
lated according to: 



Pi = TvR^r^cUl 



(2) 



where U' is the energy density of the ith component in the comov- 
ing frame. 



4.1 Constraints on the accretion luminosity and black hole 

mass 

For calculating the luminosity of the broad lines, we have followed 
Celotti, Padovani & Ghisellini (1997), namely we have assumed 
that if we set the Lya line contribution equal to 100, the total Lblh 
is 555.76, and the relative weight of the Ha, H/3, Mgll and CIV 
lines is 77, 22, 34 and 63, respectively (e.g. Francis et al. 1991). 
The information found are summarized in Tab. |4l reporting also, 
when available, the estimate of the black hole mass. When only one 
emission line is seen (as in the majority of cases, see Tab.|4) the es- 
timate of the entire BLR luminosity is uncertain. Furthermore, the 
detection of the most prominent line, the Lya one, for relatively 
nearby objects is not possible from the ground, and requires ul- 
traviolet observations from space. Plan, Falomo & Treves (2005) 
have studied a small sample of blazars spectroscopically observed 
with the Space Telescope, and compared the relative strength of the 
UV lines with the compilation of Francis et al. (1991). They found 
that the weights of Mgll and CIV are 19 and 53 (setting the line 
Lya=100), somewhat less than in Francis et al. (1991). Therefore 
the estimates given here for our blazars are uncertain by at least a 
factor 2. Despite this uncertainty, the knowledge of the BLR lumi- 
nosity gives an important constrain to the model, since it indicates 
the luminosity of the accretion disk, that we set to Ld ~ IOLblr- 
This is especially valuable when we do not have any sign of thermal 
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65 


0808+019 


2007-2009^^ 
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0.65/- 


30 
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57.6 


5.47 
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33.6±0.5 


1.0/- 


112 


1012+0630 


20109. ~ 


4.8 


1.97 


2.0±0.3 


4.6±0.6 


-/47 


45 


1026-1748 


24-07-2010 


0.01 


6.42 
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< 30.^ 


-/- 


- 


1204-0710 


09-08-2010 


5.1 


2.02 


2.8±0.3 


10.0±'o"9 


0.99/- 


4 


1338+40 


2008-2009.^ 


11.1 


0.822 


1.88±0.04 


89±2 


1.3/- 


65 


1519-273 


20-01-2010 


2.2 


9.11 


1.4±0.8 


4.2±1.2 


-/19 


11 


1807+698 


2007-2009i 


36.3 


4.11 


1.84±0.05 


31.8±0.6 


0.96/- 


76 


2007+77 


14-06-2009" 


6.3 


8.39 


1.4±0.2 


22±1 


1.08/- 


9 


2214+24B 


201o£ 


8.0 


5.75 


1.9±0.2 


10.7±0.8 


0.67/- 


7 


2240-260 


2008-2009^ 


5.5 


1.35 


1.9±0.3 


5.8±0.8 


-/55 


66 


2340+8015 


2009_;^ 


10.8 


14.2 


2.5±0.2 


9.0±0.6 


0.76/- 


9 



Table 2. Summary of XRT observations. The observation date column indicates the date of a single snapshot or the years 
during which multiple snapshots were performed. The corresponding note reports the complete set of observations integrated. 
The column "Exp" indicates the effective exposure in ks, while N^i is the Galactic absorption column in units of [10^" cm~^] 
from Kalberla et al. (2005). F is the photon index of the power law model [F{E) oc E~^}, ^02— lokoV observed 
(absorbed) flux. The two last columns indicate the results of the statistical analysis: the last column contains the degrees of 
freedom, while the last but one column displays the reduced or the value of the likelihood (Cash 1979). 



Data from different observations were integrated: 26-05-2005, 02-02-2008, 07-02-2008, 08-02-2008, 13-02-2008, 05-03- 
2010, 08-03-2010, 19-06-2010, 23-06-2010, 05-07-2010, 09-07-2010, 13-07-2010. 
^ Data from different observations were integrated: 18-05-2007, 27-02-2010. 

^ Data from different observations were integrated: 20-12-2007, 23-12-2007, 14-09-2008, 19-09-2009. 

Data from different observations were integrated: 23-10-2006, 12-12-2007, 18-09-2009,20-09-2009, 10-12-2009, 13-12- 
2009, 11-01-2010, 08-02-2010. 

^ Data from different observations were integrated: 04-07-2006, 28-03-2007, 10-01-2008, 11-01-2008, 15-01-2008, 09-01- 
2009, 01-11-2009, 05-11-2009, 12-12-2009, 23-01-2010, 05-03-2010, 12-03-2010. 

/ Best fit with a broken-power law model (ftest > 99.99%): Ti = 1.1 ± 0.1, Ebrcak = 1.3 ± 0.1 keV, r2 = 1.89±0.07. 
In the table is indicated r2 only. 

9 Data from different observations were integrated: 24-05-2010 (two observations), 25-05-2010. 

^ Upper limit derived with PIMMS with fixed photon index equal to 2. 

* Data from different observations were integrated: 15-10-2008, 21-12-2009. 

Data from different observations were integrated: 01-03-2007, 15-04-2007, 22-01-2009, 09-1 1-2009, 1 1-1 1-2009. 
^ Data from different observations were integrated: 25-01-2010, 27-01-2010. 
' Data from different observations were integrated: 25-12-2008, 22-09-2009. 

Data from different observations were integrated: 05-09-2009, 08-09-2009, 09-09-2009. 



emission in the opdcal-UV, often dominated by the non-thermal 
continuum. In these cases we have also chosen a value for the black 
hole mass consistent with what found in the literature. 



4.2 Results of the modelling 

In Fig.[T]-[8]of the Appendix we show the SED of the considered 
BL Lacs and the fitting model. The parameters for the modelling 
are listed in Tab. [5] and the derived jet powers in Tab. 6. One key 
property of the majority of our sources is to have a relatively low lu- 
minosity accretion disk. If the size of the BLR is connected with Ld 
1/2 

(we assume i?BLR oc ) then this implies very compact sizes of 
the BLR, both in absolute terms and in units of the Schwarzschild 
radius. On the contrary, the dissipation region is always at a few 
or several hundreds of Schwarzschild radii, and in 24/30 cases we 
have i?diss > i^BLR. This is in agreement with what found in GIO, 
but here this issue can be treated in more detail thanks to the knowl- 



edge, for some of the sources, of the luminosity of some broad lines 
and the black hole mass. In Fig. [T]-[8]of the Appendix we show, 
separately, the contribution of the synchrotron self-Compton (SSC, 
long dashed line) and of the External Compton (EC, dot~dashed 
line) components. We find a variety of cases, from sources whose 
high energy bump is completely dominated by the EC component 
(see e.g. 0058+331 1; 0208-512; 1803+784), or by a pure SSC (e.g. 
0521-365; 0558-3838; 0851+202; 0907+3341; 1026-174; 1057- 
79; 1147+24; 1204-71; 1557+565; 2340+801), or by the SSC at 
softer X~ray energies and by the EC at higher energies. Rarely (see 
0954+658 and 2200+420) there is an important contribution by the 
second order Compton scattering in the SSC spectrum, competing 
with the EC component in the GeV band. We alert the reader that 
in some cases (for instance: 0558-3838 and 1028-174) the paucity 
of the data points makes the resulting "fit" very uncertain. 

There are 8 sources in common with the sample studied by 
GIO (i.e. all blazars with redshift in the LBAS catalog). These 
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Table 3. Summary of Swift/UVOT observed magnitudes. Lower limits are at 3cr level. 



sources are indicated in italics both in Tab.[T]and in Tab.|5] With re- 
spect to the model parameters adopted in GIO, we here have talcen 
into account the luminosity of the BLR and the existing estimates 
of the black hole mass. 

For 0208-512, the black hole mass adopted here is AI — 
10^ Mq (it was 7 X 10* M© in GIO), T = 13 (it was T = 10), 
P' is nearly half than in GIO, Ld is similar, i? = 3 G (it was 2.05). 

Also for 0537^41 we have changed the adopted black hole 
mass (now M = 6 x 10** Mq, it was 2 x IO^Mq in GIO), halving 
the accretion disk luminosity. Note that these changes in black hole 
mass, that are rather large, are possible in sources whose IR-UV 
continuum is dominated by the non-thermal beamed emission. In- 
stead, when the accretion disk emission dominates, the estimate of 
the black hole mass is much less uncertain (see the discussion in 
GIO and Ghiselhni et al. 2010b). 

For 0426-380 the estimated BLR luminosity implied a much 
reduced accretion disk luminosity with respect to what assumed in 
GlO, and this in turn implied a smaller i?diss, larger magnetic field, 
and larger P. Instead, for 1057-79, the smaller value for resulted 
in -Rdiss slightly larger than in GlO and in a smaller magnetic field. 

For 0851+202 (=0J 287), 1749+096 and 2200+420 (=BL Lac) 
we can now fix the value of Ld (only an upper limit was used in 
GlO, consistent with the values used now). For these sources the 
model parameters are quite similar to the ones used in GIO. 

For the sources analyzed in this paper and not present in LBAS 
the parameters are similar to the ones derived in GIO for the blazars 
of similar 7-ray luminosity. The location of the dissipation region 
7?diss is a several hundreds Schwarzschild radii, the bulk Lorentz 
factor is in the range 10-15 (with the exception of 0521-365, with 
r — 5), magnetic field B ~0.25-7 G, black hole masses M ~(3- 
10) X 10* Mq . The possible, expected, difference with respect to the 
blazars in GIO is that sometimes the viewing angle 6^ is larger than 
the typical value of 3 degrees (see 0521-365 with 9^ = 12°). This 
is expected because, with the smaller limiting flux of the ILAC 
catalog (with respect to LBAS) we start to see not only the sources 
maximally Doppler enhanced (i.e. pointing almost exactly at us), 
but also the blazars that are slightly misaligned. 

The main point that the modelling should help to clarify is why 
these intruder BL Lacs have a relatively steep spectrum, more sim- 
ilar to the slopes of FSRQs than to the rest of BL Lacs. To answer 
this issue, consider first those intruder BL Lacs that are instead FS- 



RQs, with luminous accretion disk, dissipation regions within the 
BLR, and therefore characterized by 7-ray slopes similar to "pure" 
FSRQs. In these sources radiative cooling is very severe (see the 
last column of Tab.|5]reporting the values of 7c), and the 7-ray flux 
is always produced by the steep part of the electron distribution. 

The remaining objects, instead, have low luminosity accretion 
disks, small -Rblr, and dissipation locations beyond the BLR (i.e. 
-Rdiss > ^^blr)- Nevertheless, contrary to BL Lacs showing a flat 
spectrum (i.e. a-y < 1), they also are characterized by a relatively 
strong radiative cooling, in this case mainly due to the synchrotron 
and the SSC processes. As a result the electron population produc- 
ing the 7-ray flux is steep (i.e. a-, ~ S2/2 above 7c). Finally, in 
a few cases (see for instance 1807+698=3C 371), the cooling is 
weak, with 7c becoming equal to 7max, and the 7-ray spectrum is 
steep because it is produced by most energetic electrons (near the 
end of the electron energy distribution). 

We would like to stress that, apart from finding the reason of 
the steep 7-ray spectrum, the aim of the present paper is not on the 
outcomes of the model fitting, being instead the finding of a new 
classifying scheme for the sub-classes of blazars. However, know- 
ing the intrinsic physical quantities characterizing the emitting re- 
gion of the jet is certainly helpful for our goal, as is the determina- 
tion of the black hole mass and accretion rate of those blazars with 
the IR-UV continuum dominated by the accretion luminosity. 



5 A NEW CLASSIFICATION FOR FLAT SPECTRUM 
RADIO QUASARS AND BL LAC OBJECTS 

Tab. |4] lists the luminosities and EW of the broad lines of our in- 
truder BL Lacs. For several of these objects there are estimates of 
the mass of their black holes and therefore we have the luminos- 
ity of the BLR in Eddington units. For the sources with no black 
hole mass estimate, we use a fiducial, average, black hole mass 
Af = 3 X 10* Mq. 

The last column of Tab. |4] reports a tentative classification of 
the object on the base of its SED, independent of the presence of 
the broad emission lines and of their strength. For this SED-based 
classification we follow following criteria: 

• "FS" are those sources, currently classified as BL Lacs, whose 
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Table 4. Emission lines, BLR total luminosities, black hole masses, and BLR luminosities in units of Eddington ones. Col [3] reports the maximum observed 
equivalent width EW in A. In Col. [8] the number in parenthesis is the value of the black hole mass used. When the black hole mass is unknown, we have 
assumed logMeH/JV/Q = 8.5. For 0235+164 and 0426-380 we have used logMeH/A/© = 9 in agreement with our previous estimates (Ghisellini 
et al. 2009) derived from fitting the SED. The last column gives the classification according to the appearance of the SED shown in Figs. |2]tl] and of the 
presence/absence of prominent broad lines. Question marks mean that the classification is uncertain. 17/30 (57%) have broad lines; 6/30 (20%) are "pure" 
FS; 17/30 (57%) are LBL; 3/30 (10%) are IBL; 4/30 (13%) are HBL. The last three entries are BL Lacs present in GIO for which we found data for the 
broad emission lines. All these thi'ee are HBL. References for emission lines: Ca03: Carangelo et al. 2003; CG97: Celotti, Padovani & Ghisellini 1997; La96: 
Lawrence et al. 1996; LaOl: Landt et al. 2001; Pi05: Plan, Falomo & Treves 2005; Ra07: Raiteri et al. 2007; ReOl: Rector & Stocke 2001; Sc97; Scarpa & 
Falomo 1997; Sb05: Sbarufatti et al. 2005; Sb06: Sbarufatti et al. 2006; Sb06: Sbarufatti et al. 2009; SDSS: http://cas.sdss.org St89: Stickel et al. 1989; St93: 
Stickel et al. 1993; Ve95: Vermeulen et al. 1995; White et al. 1988. References for the black hole masses: Ba03: Barth et al. 2003; Fa03: Falomo et al. 2003a; 
Fa03b: Falomo, Carangelo & Treves 2003b; Fa04: Fan & Cao 2004; Li06:Liu et al. 2006; Wa04: Wang, Luo & Ho 2004; Wa08: Wagner 2008; Wo05: Woo & 
Urry 2005; Wu02: Wu, Liu & Zhang 2002. 



high energy peak dominates the bolometric output and the X-ray 
spectrum belongs to the high energy peak; 

• "LBL" (i.e. Low peak BL Lacs, according to the definition 
of Padovani & Giommi 1995) have the synchrotron peak and high 
energy peak of comparable luminosities, and the X-ray emission 
belongs to the high energy peak; 

• "IBL" (i.e. Intermediate BL Lacs) are defined according to the 
shape of the X-ray spectrum, steep (i.e. due to the tail of the syn- 
chrotron emission) at low frequencies and hardening at higher fre- 
quencies; 



• "HBL" (High peak BL Lacs) have the X-ray emission dom- 
inated by the synchrotron process. These objects have usually a 
steep X-ray spectrum, but sometimes, as in the famous case of the 
flaring state of IVIkn 501 (Plan et al. 1998), the synchrotron spec- 
trum peaks at so large frequencies to make the X-ray spectrum flat. 

The above distinctions builds on the classification proposed by 
Padovani & Giommi (1995; 1996) (see also the recent extension 
in Abdo et al. 2010b) on the base of the position of the synchrotron 
peak energy, and it is consistent with the typical SEDs we have 
studied in GIG and in Tavecchio et al. (2010). Of the 30 objects 
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Table 5. List of parameters used to construct the theoretical SED. Not all of them are "input parameters" for the model, because /Jblr is uniquely determined 
from Ld, and the cooling energy 7c is a derived parameter. Col. [1]: name; Col. [2]: redshift; Col. [3]: dissipation radius in units of 10^^ cm and (in parenthesis) 
in units of Schwarzschild radii; Col. [4]: black hole mass in solar masses; Col. [5]: size of the BLR in units of 10^^ cm; Col. [6]: power injected in the blob 
calculated in the comoving frame, in units of 10*^ erg s~^; Col. [7]; accretion disk luminosity in units of 10^''' erg s^^ and (in parenthesis) in units of LEdd; 
Col. [8]: magnetic field in Gauss; Col. [9]: bulk Lorentz factor at Rd{ss\ Col. [10]: viewing angle in degrees; Col. [1 1] and [13]: minimum, break and maximum 
random Lorentz factors of the injected electrons; Col. [14]: and [15]: slopes of the injected electron distribution [(3(7)] below and above 7b; Col. [16]: values 
of the minimum random Lorentz factor of those electrons cooling in one light crossing time. The total X-ray corona luminosity is assumed to be in the range 
10-30 per cent of Lj. Its spectral shape is assumed to be always oc exp(— /ii//150 l;eV). 



that we considered, 17/30 (57%) have detected broad Unes; 6/30 
(20%) are FS; 17/30 (57%) are LBL; 3/30 (10%) are IBL; 4/30 
(13%)areHBL. 

The relation between I/BLR/iEdd and the 7-ray luminosity 
(measured in Eddington units as well) is shown in Fig. [3] where 
we plot the objects with detected broad lines, divided according 
to the classification from the SED. To these sources we have also 
added all the blazars detected in the 3-months all sky survey of 
FermifLKT (LB AS, see A09 and GIO) for which there are estimates 
of the mass of their black holes and flux measurements of at least 
one broad emission line. These are 3 HBL (i.e. Mkn 421, Mkn 501 
and 2005^89; data in Tab.© and 18 FSRQs (data in Tab. IT}. 

A clear trend is visible: objects with stronger emission lines 
are more luminous in the 7-ray band (normalizing both luminosi- 
ties to the Eddington one). In this figure we have drawn a dividing 
line at Lblr/ L^dd = 5 x 10^^, that separates "pure FSRQs" and 
"FS" from the rest. The considered sources are still limited in num- 
ber, but this division is in agreement with the idea that the blazars' 
divide occur for a change in the accretion regime. If we consider 



a disk luminosity roughly 10 times greater than Lblr, the separa- 
tion from pure FSRQs to BL Lacs occurs at Ld ~ 5 x 10~'^LEdd, 
approximately where the disk goes from a radiatively efficient to 
inefficient regime. This is in remarkable good agreement with what 
we suggested earlier (GMT09 for Fermi blazars and Ghisellini & 
Celotti 2001 for FRI and FRII radio-galaxies). 

BL Lac objects that we have reclassified as "FS" indeed oc- 
cupy the same region of FSRQs, confirming that the EW-based 
classification scheme sometimes hides the real nature of the source. 
"LBLs" are intermediate sources, where the line luminosity (in Ed- 
dington units) decreases (but it does not disappear) as well as the 
7-ray luminosity. "HBLs" are at the extreme of the distribution, 
with very weak lines and weak L-y. Consider furthermore that Fig. 
[3] shows only the sources with measured broad lines (there is only 
one exception, 1204-071, for which we have found an upper limit), 
and therefore we should not give too much weight to the apparent 
correlation between I/BLR/^Edd and L^/L-Edd- Consider also that 
luminosity variations exceeding even two orders of magnitude are 
not uncommon for these objects. Bearing these caveats in mind, 
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Table 6. Logarithm of the jet power in the form of radiation (Pi), Poynting 
flux (Pb). bulk motion of electrons (Po) and protons (Pp, assuming one 
proton per emitting electron). Powers are in erg s^^. 



we nevertheless believe that this figure indeed suggests that the se- 
quence HBL IBL (hopefully, no IBL is included, yet) — > LBL 
— >■ FSRQ can be explained as a sequence of strength of the broad 
lines. 

The location of the points seems to suggest a continuity of 
isLR/iEdd values, rather than a bimodal distribution, but the 
number of sources is really too small for a strong claim. We then 
leave this issue to future studies: these will give some insight on the 
properties of the ionizing flux: is it changing dramatically when the 
disk emits less than a few thousandth of the Eddington luminosity, 
or smoothly? 

An important caveat is in order: we considered only the 
sources of known redshift, and as discussed in §2 there are many 
sources classified as BL Lacs in the ILAC sample of unknown red- 
shift. In order to discuss how these sources can affect the proposed 
classification scheme, consider the possibility that they all are at 
z — 0.5 or else at z = 2. In the first case, as Fig. [T] shows, they 
would have a low or moderate 7-ray luminosity. The fact that their 
spectrum is featureless implies that their emission lines are intrin- 
sically weak, and if their black hole mass is similar to the other BL 
Lacs their LblhJ L^dd and L^/L^dd ratios locate these objects in 
the bottom left region of Fig. [3] together with the shown sources. 

Instead, if z ~ 2, these sources would have a 7-ray luminosity 
as large as the powerful FSRQs (the subset of these sources having 
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Table 7. The luminosity of the BLR in units of Eddington for the 7-ray 
FSRQs studied in GIO, for which there is an estimate of the black hole mass 
M. The specific value used for log M is indicated in square brackets in the 
last column. " : the SDSS associates the Mgll line with the Alll+Clll line, so 
the SDSS redshift is wrong. References: Be02: Bechtold et al. 2002; Ce97: 
Celotti, Padovani & Ghisellini 1997; Li06: Liu et al. 2006; Os94: Osmer et 
al. 1994; Pi05: Plan, Falomo & Treves 2005; SDSS; Schneider et al. 2010, 
Shen et al. 2010 (https;//www.cfa.harvard.edu/~yshen/BH_mass/dr7.htm ); 
St93: Stickel et al. 1993. 



•o 
u 



m 
J 

OO 
O 



-2 - 



-4 - 



~i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



o 



O FSRQs 

• 'FS' 

* LBL 
□ HBL 



O 



O 



O 



o o 

QD 

* 



•k 



□ 

I I I I 



■k 



-3-2-10 1 2 

Log LyLj,dd 

Figure 3. Luminosity of the broad line region (in units of the Eddington 
one for sources with at least one broad line in their spectrum and with an 
estimate of the black hole mass) as a function of the 7-ray luminosity in 
units of the Eddington one. Empty (red) circles: FSRQs studied in GIO 
with estimates of the black hole mass and broad emission line luminosity; 
full (violet) circles: TS' sources (i.e. BL Lacs reclassified as FSRQs in this 
paper); (green) stars: LBLs; empty (blue) squares and upper limit: HBLs. 
The three data points for HBLs (Mkn 421, Mkn 510 and 2005^89) do not 
belong to the sample studied in this paper (i.e. they have a-y < 1/2), but 
are included for comparison. 
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> 1.2 indeed occupies the same region occupied by FSRQs in 
Fig.[TJ. We iiave two possibilities: 

• They have powerful lines and an even more powerful contin- 
uum, swamping the emission lines. In this case, and for average 
black hole masses, their I/BLR/iEdd is large, as it is L-y/L^dd- 
They would then occupy the top right part of Fig. [5] together with 
the other FSRQs. They would be sources similar to 0208-512, but 
with an even stronger continuum. 

• They have intrinsically weak lines, and average black hole 
masses. In this case I/BLa/iEdd would be small, but L^/L-^dd 
is large. They would then occupy the bottom right part of Fig. [3] 
now devoid of sources. In this case there would be no correlation 
between LBLR/^Edd and L-y/L-Edd- These sources would have a 
powerful non-thermal emission, possibly extremely beamed (sug- 
gesting a very powerful jet), even if their accretion disk is weak. 

We conclude that the new classification scheme we are propos- 
ing is not affected by the possible redshift values of the BL Lacs 
with unknown z. On the contrary, our physically-based scheme 
could strengthen possible problems concerning the blazar se- 
quence. In fact, in the latter mentioned case (i.e. z ~ 2 for the 
sources of unknown z) we could demonstrate the existence of very 
powerful sources that are classified as BL Lacs not only according 
to the "old" definition (based on the EW), but also according to the 
one we are proposing. 

The suggested new classification scheme can have important 
consequences and ramifications, that we discuss below. 



6 DISCUSSION 

The most important result of this study is the suggestion of a 
new classification criterion distinguishing BL Lac objects from FS- 
RQs, based on the luminosity of the broad emission lines mea- 
sured in Eddington units. The critical value we propose is around 
iBLR/i/Edd ~ 5 X 10"*. We base this suggestion on the mea- 
surements of the emission line luminosities, on the estimates of 
the black hole mass and on the conversion of the luminosity of a 
few lines (often, only one line) into the luminosity of the ensem- 
ble of the broad emission lines. We find that the I/BLR /iEdd ratio 
is proportional (with some scatter) to the 7-ray luminosity mea- 
sured in Eddington units. If we furthermore divide the blazars into 
sub-categories according to the properties of their SED, then again 
we find a trend: low values of the I/BLR/ifidd ratio correspond to 
low power, high energy emitting (i.e. "blue") BL Lacs. Vice-versa, 
high values of I/BLR/I^Edd correspond to powerful "red" FSRQs. 
This is in total agreement with the blazar sequence, and also with 
the blazars' divide we have proposed earlier. 

We would like to stress that the distinction between FSRQs 
and BL Lacs we are proposing here does not imply a dichotomy. 
This could be present if the entire accretion disk suddendly changes 
regime below and above a given threshold in M, but more likely 
there should be a gradual change (namely, some parts of the accre- 
tion disks might be radiatively efficient and others not). Then the 
transition between BL Lacs and FSRQs might be gradual as well. 

The first example of how this new classification might be use- 
ful concerns our starting sample of "intruder" BL Lacs. Fig. |4] 
shows how the different types of intruder BL Lacs are located in 
the a-y-L^, plane. Blazars re-classified as "FS" all have large 7- 
ray luminosities (filled magenta circles) and fall in the same region 
of luminous FSRQs. LBLs span a large range of L^, but not the 
very high end, and can be thought as intermediate objects. HBLs 
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Figure 4. Same as Fig.[T] but with the intruder BL Lacs reclassified accord- 
ing to Tab.|4] Filled (violet) circles: sources reclassified as FSRQs; (green) 
stars: LBLs; (cyan) filled squares: IBLs; (blue) diamonds: HBLs. 



are at significantly lower 7-ray luminosities. Of course one is left 
curious to see what happens for the remaining BL Lacs of this plot. 
To satisfy this curiosity we plan to systematically study them (Sbar- 
rato et al., in preparation). 

The other obvious advantage of our proposed classification is 
that it is physically based, and will help to construct cleaner sam- 
ples of objects aiming to study possible different properties be- 
tween BL Lacs and FSRQs, or, rather, between weak and strong 
line objects. For instance, this might help to clarify why some BL 
Lac objects, observed and imaged in the radio band, appear to be 
FR II radio-galaxies (see e.g. Kharb, Lister & Cooper 2010) con- 
trary to the accepted scenario of the parent population of blazars, 
assigning FR I to BL Lacs and FR II to FSRQs. 

A particularly interesting issue that the new classification 
scheme will help to clarify is the evolution of blazars. Hints of 
no or slightly negative evolution of BL Lacs (or subcategories of 
them) and of positive evolution of FSRQs could be associated with 
the evolution of the accretion rate in cosmic time. It is possible 
that a better understanding may come from considering the entire 
blazar population as a whole, characterized by larger rates of accre- 
tion in the past (and therefore by a prevalence of FSRQs over BL 
Lacs) and a decreased rate of accretion now (with BL Lacs becom- 
ing more numerous). We believe that this idea, put forward, among 
others, by Maraschi & Rovetti (1994) and by Cavaliere & D'Elia 
(2002) is worth pursuing. 

By linking the BLR luminosity to the underlying accretion lu- 
minosity we suggest that the different "look" of blazars reflects 
primarily their accretion rate in units of the Eddington one. The 
M / M-Edd ratio controls the radiative efficiency of the accretion 
disk, and thus the absolute luminosity of the broad lines, and 
(through the BLR size - disk luminosity relation), their distance 
from the black hole. This in turn regulates the importance of the 
External Compton process for the formation of the high energy 
bump of the SED. Powerful blazars do have powerful disks, lumi- 
nous broad lines emitted at large distances. If jet dissipation occurs 
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Figure 5. Power of the jet spent in the for of radiation P,- as a function of the 
accretion luminosity . Black symbols are the estimates in this paper, BAT 
points (grey diamonds) come from the high-redshift blazars present in the 3 
year all sky survey of BAT (and studied in Ghisellini et al. 2010b), the other 
points and upper limits come from GIG, and are divided according to their 
7-ray luminosities, as labelled. The grey stripe indicates what expected if 
Ljj oc Al for all luminosities, while oc and oc Al at low and 
high values of M, respectively. 



ones studied in GIO and in Ghisellini et al. (2010b). For large values 
of Ld and Pr, the two quantities are proportional. Instead, below 
Ld ~ 10*^ erg s^^, the data points and the shown upper limits are 

1/2 

consistent with Pr oc . 

We interpret this behavior (see also GIO and Ghisellini & 
Tavecchio 2008) as follows. 

Assume that Pr is proportional the total jet power Pjct, in turn 
always proportional to the accretion rate Mc^. Instead, assume that 
Ld ~ rfMc? , with rj being constant only above a critical luminos- 
ity (in Eddington units), while below this critical value rj oc M 
(Narayan, Garcia & McClintock 1997). This means that M oc Ld 



above some critical value Lc, and A4 oc 
Pr oc Al thus implies 



1/2 



Pr 
Pr 



OC 
OC 



-1/2 



Ld > Lc 
Ld < Lc 



below. Assuming 



(3) 



If the jet power is always proportional to M, its intrinsic prop- 
erties should not change according to Al being smaller or greater 
than any critical value. In other words the jet properties should not 
depend on M/MEdd- Therefore the jet power traces Al better than 
the accretion luminosity, that strongly depends upon M /M-Edd ■ On 
the other hand, the "look" of the jet (i.e. the produced non-thermal 
SED) is strongly influenced by A//A/Edd because it is the latter 
ratio that rules the strength of the external radiation used as seed 
for the formation of the high energy bump of the SED. 

To reiterate: the jet power does depend on Al linearly, but 
there is no dependence of the formation, collimation and acceler- 
ation mecahnisms of the relativistic jet in blazars on A//MBdd- 
Relativistic jets are present for all values of A//A/Edd- 



within the BLR, the External Compton process is dominant and 
the cooling severe, resulting in 7-ray dominated sources. Weaker 
blazars have weak disk and weaker lines emitted closer to the black 
hole. Dissipation in the jet occurs outside the BLR, the main pro- 
cesses becomes synchrotron and SSC, resulting in more equally 
shared luminosity between the synchrotron and the SSC compo- 
nents. 

The other strong evidence we have found in our earlier works 
(e.g. Celotti & Ghisellini 2007; GIO; Ghisellini et al. 2010b), is 
that the total jet power Pjet is proportional to Al, and indeed very 
close to Mc^ independently of the accretion regime. The claim that 
Pjet ~ Mc^ depends somewhat to the assumption to have one 
proton per emitting electron. Although we do have limits to the 
possible amount of pairs (that cannot be energetically important, 
see Sikora & Madejski 2000, Celotti & Ghisellini 2008; Ghisellini 
& Tavecchio 2010) some uncertainty remains on the exact amount 
of pairs in jets. An almost model-independent lower limit on the 
jet power is the power spent by the jet to produce the non-thermal 
radiation, Pr. It is a lower limit because if Pr were comparable to 
the total jet power, Pjot, then the jet would strongly decelerate, and 
no superluminal motion would be seen. Assuming that Pr is pro- 
portional to Pjct implies a constant efficiency in converting part of 
the bulk relativistic motion into relativistic electrons and then ra- 
diation. Pr can be derived by Eq. 2 substituting, for (7i, the energy 
density of the radiation produced by the jet, evaluated in the comov- 
ing frame. The only free parameter entering in the estimate of Pr 
is the bulk Lorentz factor (Pr oc L/F^, where L is the bolometric 
luminosity of the jet, calculated assuming isotropy). 

In Fig.[5]we show Pr as a function of the accretion luminosity 
Ld- It shows the blazars analyzed in this paper together with the 



7 SUMMARY AND CONCLUSIONS 

In this work we have analyzed a sample of blazars detected by 
FermifLAT that have been classified as BL Lac objects and that 
have an energy 7-ray spectral index a.y > 1.2. They therefore oc- 
cupy a region, in the ot~,-L-y plane, preferentially occupied by FS- 
RQs. Our intent was to investigate the properties of these objects, 
to see if they could be considered as intermediate objects between 
"pure" BL Lacs and "pure" FSRQs. 

Doing so, we collected from the literature the broad emission 
line data for a sizeable number of these sources, as well as estimates 
of their black hole mass. At the same time, we model their SED 
with a one-zone leptonic model, to find out their intrinsic properties 
and especially to investigate why their 7-ray spectrum is rather 
steep. Our main results are the following: 

• From the model fitting, we explain the relatively steep a-y of 
these blazars as due to a relatively severe cooling of the electron 
population. The cooling is particularly severe in sources that have 
strong disks and emission lines, but also in the remaining sources 
it is fast enough to make the emitting electron distribution steep. 

• Some of the considered blazars, classified as BL Lacs, have 
broad emission lines as strong as in FSRQs, both in absolute terms 
and in Eddington units. 

• There is a trend associating the BLR luminosity in Eddington 
units with the 7-ray luminosity. Due to the paucity of points, we 
cannot claim that there is a strict correlation, yet the indication is 
that I/BLR/LEdd, I/7/I/Edd and the type of the SED (i.e. LBL or 
HBL) are strongly linked. 

• From this evidence, we suggest a new classification scheme 
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for dividing BL Lacs from FSRQs, based on the BLR luminos- 
ity in Eddington units: we propose to set the dividing value at 

iBLR/^Edd ~ 5 X 10"*. 

• Since the BLR is thought to intercept and reprocess about 
10% of the disk luminosity, the dividing value corresponds to a 
disk emitting at the ~0.5% of the Eddington limit. This is, approx- 
imately, also the value dividing the radiatively efficient from the 
radiatively inefficient accretion regimes. 

• This work, together with the previous studies we have done 
on Fermi (and EGRET) blazars, confirms that jets are powerful, 
and that they are born and laimched for all values of the accretion 
rate (in Eddington units). 
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APPENDIX: SPECTRAL ENERGY DISTRIBUTION 



© 0000 RAS, MNRAS 000, 000-000 



14 



G. Ghisellini et al. 



■ 1 1 



15 



_ 1 1 1 1 1 1 1 

_ 005B+3311 


1 1 1 1 1 1 




1 _ 
- 


\ z=1.371 








- 




EC / 




- 




o 










disk 

i 












li 






- 8 1 

- , 1 1 


tpriis^ \ = 

\ 1 X 

J , /, iT 


/ 
1 

1 corona 

r--,---,.. , 1 , 


^\ 


1 - 


10 


15 


20 




25 



47 




15 20 
Log u [Hz] 

Figure 1. SED of 0058+331 1 and 0109+22. Darker points (red in the elec- 
tronic version) refer to tlie fermi/LAT and Swift (UVOT and XRT) obser- 
vations. The lines are the result of the modelling. We label the synchrotron 
component (green solid Une), the disk, torus and X-ray corona emission 
(dotted black), the SSC flux (grey long dashed) and external Compton flux 
(grey dot-dashed). The thick (blue) line is the sum of all components. 
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Figure 2. SED of PKS 0208-5 12, PKS 0235+164, PKS 0426-380 and PKS 
0521-365. Symbols and lines as in Fig.|2] 
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Figure 3. SED of PKS 0537^41, PMN 0558-3839, PKS 0754+100 and Figure 4. SED of PKS 0829+046, 0851+202 (=OJ 287), 0907+3341 

PKS 0808+019. Symbols and lines as in Fig.|2] (=TON 1015) and 0954+658. Symbols and lines as in Fig.|2] 
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